Discrete antiferromagnetic spin-wave excitations in the giant ferric wheel Fei 8 
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The low-temperature elementary spin excitations in the AFM molecular wheel Fei8 were studied 
experimentally by inelastic neutron scattering and theoretically by modern numerical methods, 
such as dynamical density matrix renormalization group or quantum Monte Carlo techniques, and 
analytical spin-wave theory calculations. Feig involves eighteen spin-5/2 Fe m ions with a Hilbert 



space dimension of 



10 1 



constituting a physical system that is situated in a region between 



microscopic and macroscopic. The combined experimental and theoretical approach allowed us to 
characterize and discuss the magnetic properties of Feis in great detail. It is demonstrated that 
physical concepts such as the rotational-band or L&i?-band concepts developed for smaller rings 
are still applicable. In particular, the higher-lying low-temperature elementary spin excitations in 
Feis or AFM wheels in general are of discrete antiferromagnetic spin-wave character. 



PACS numbers: 75.50.Xx, 75.10.Jm, 78.70.Nx 



I. INTRODUCTION 

Ring-like arrangements of a dozen or so of magnetic 
spins experiencing nearest-neighbor antiferromagnetic 
(AFM) exchange interactions, as realized experimentally 
for instance by the AFM molecular wheels, have at- 
tracted significant attention in the past decade The 
molecular ferric wheel [Fei8(pdH)i2(0 2 CEt) 6 (N03) 6 ], or 
Feis hi short, is the largest magnetic molecular wheel 
synthesized to date^ The molecule contains N = 18 
Fe m ions with spin s = 5/2, arranged in a ring- like fash- 
ion as shown in Fig. [IJa) , and exhibits crystallographic 
Cg symmetry. Its large yet finite size makes it an ideal 
candidate to explore the region between microscopic and 
macroscopic physics (the system is mesoscopic). The 
magnetism in the Feis wheel was studied before using 
high-field magnetic torque measurements, and the low- 
lying energy spectrum up to 2 meV was investigated by 
inelastic neutron scattering (INS)J^ The experimental 
data demonstrated the dynamics of the Neel vector, and 
the magnetic torque provided direct evidence for quan- 
tum oscillations in the Neel vector tunneling gap due to 
quantum phase interference^ 

The experimental observations have been well de- 
scribed in terms of an effective two-sublattice Hamilto- 
nian, which had been demonstrated before for smaller 
wheels with up to 10 spin sites to approximate well the 
low-energy part of the true spectrum ! 15 : 16 This Hamil- 
tonian is in fact related to the more general concept 
of the L&i?-band picture, which describes the elemen- 
tary excitations as a set of rotational (parabolic) energy 
bands, and which was shown to apply to a variety of 
AFM spin clusters with bipartite or tripartite sublat- 
tice structure i 10 ' 1 5 ~ — The effective Hamiltonian allowed 



a description of the experiments because it operates in 
a Hilbert space of dimension 2116 (for Feis), which can 
easily be handled on a personal computer using stan- 
dard numerical diagonalization techniques. However, di- 
rect confirmation of its applicability to wheels as large as 
Feis is lacking and important magnetic parameters such 
as magnetic anisotropy could not be determined reliably. 
Furthermore, higher-lying spin excitations, which are ex- 
pected in the L&E'-band concept, were not observed. 

The numerically exact evaluation of all energy eigen- 
values and eigenfunctions of a giant molecule such as 
Feis poses a great challenge for theory since the size 
of the related Hilbert space grows as (2s + 1)^ and for 
Feis assumes a value of 101,559,956,668,416 10 14 . 
This dimension is much too big for a matrix diagonal- 
ization. Even a decomposition of the Hamiltonian ma- 
trix according to the available symmetries, as success- 
fully done for the smaller AFM wheels such as CsFes 
or Feio^^r— is not efficient enough to ease the prob- 
lem. Fortunately, the magnetic molecule Feis is non- 
frustrated, which permits the application of Quantum 
Monte Carlo (QMC) methods^! - — and furthermore is 
quasi one-dimensional, which makes it ideal for (Dynam- 
ical) Density Matrix Renormalization Group (DDMRG 
& DMRG) calculations The first method allows the 
evaluation of thermodynamic observables such as the 
magnetic susceptibility whereas the second delivers tran- 
sition rates between low-lying energy levels that can be 
related to the INS spectrum. 

In this work a comprehensive study of the higher-lying 
excitations in Feis, which in the language of the L&zE- 
band concept correspond to the E band or discrete spin- 
wave excitations, is reported. Experimentally the ex- 
citation spectrum was determined by high-energy INS 
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FIG. 1: (Color online) (a) Ball-and-stick representation of the 
molecular structure of the Feis molecule (Fe: green, O: red, 
N: blue, C: gray, and H: white). The dashed box indicates 
the unit cell used for magnetic modeling, see next panel, (b) 
Unit cell (dashed box) and labeling of the spin sites (text at 
bottom) used for the theoretical calculations as suggested by 
the symmetry of the molecule. The exchange constants Ji, 
J2, and J3 associated to different bonds are indicated. 



measurements extending the energy range to 13.5 meV, 
and the temperature-dependent magnetic susceptibility, 
which probes the full energy spectrum. Theoretically, 
the QMC and DDMRG techniques were used to repro- 
duce the experimental magnetic susceptibility and INS 
data with excellent accuracy. Because of the size and 
the structure of Feis the DDMRG calculations are time- 
consuming, which prevents a systematic scanning of the 
magnetic parameter space or least-squares fit approaches. 
This problem was circumvented by resorting to spin- wave 
theory as an intermediate step. This allowed us to refine 
the microscopic spin Hamiltonian for Feis and deduce ac- 
curate microscopic magnetic parameters. A key result is 
that the spin- wave like character of the higher lying exci- 
tations is preserved although the refined Hamiltonian is 
less symmetric compared to the assumptions in Ref. [Til . 
Furthermore, the validity of the L&iS-band concept is 
confirmed for Feis. 
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FIG. 2: (Color online) Experimental (open circles) and sim- 
ulated (lines) magnetic susceptibility as function of tempera- 
ture. The three simulated curves were obtained from the mi- 
croscopic model Eq. fT} using QMC. The parameters in the 
simulations are as indicated (with J\ — J2 = J and J3 = J'). 

The article is organized as follows: It begins with a 
discussion of the experimental results for the susceptibil- 
ity and the INS cross section in Sec. [TTJ This is followed 
by a theoretical analysis in Sec. [TTT] and a discussion in 
Sec. IIVI After our final conclusions further details of the 
employed methods are given in an appendix. 

II. EXPERIMENTAL RESULTS 

The experimental methods used for sample prepara- 
tion, and magnetic susceptibility and INS measurements 
are described in Appendix IA II Figure [2] shows the mag- 
netic susceptibility \ as function of temperature mea- 
sured on a poly-crystalline sample. At 250 K, a x value of 
0.17 cm 3 /mol is observed, which increases with decreas- 
ing temperature, reaches a maximum of 0.23 cm 3 /mol at 
a temperature of ca. 50 K, and decreases further with 
decreasing temperature. At the lowest temperature an 
abrupt decrease is observed. Such a behavior is typi- 
cal for AFM molecular wheels^ and the data are con- 
sistent with earlier measurements on a micro-crystalline 
sample.— The ground state of an even-membered AFM 
wheel is a total spin singlet, 5 = 0, and the susceptibility 
is hence expected to approach zero at zero temperature, 
which is not observed here because in large wheels the 
drop to x = occurs only at the lowest temperatures 
(the gap to the first excited triplet is roughly given by 
~ 4J /N and is small in Feis) and the susceptibility in this 
temperature regime is strongly affected by the presence of 
magnetic anisotropy (which is significant in Feis)^ The 
solid curves are theoretical results which are discussed 
below. 

The INS spectrum recorded at an incoming wavelength 
of A = 4.2 A and a temperature of 1.9 K on the spectrom- 
eter IN5 is shown in Fig. [3] At low energies a strong 
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FIG. 3: Experimental INS spectra recorded on IN5 with in- 
coming neutron wavelength A = 4.2 A at 1.9 K. Positive en- 
ergy transfer corresponds to neutron energy loss. The labels 
indicate observed transitions. The inset sketches the assign- 
ment of the low-energy excitations as inferred previously.— 



feature at 0.3 meV (peak la) together with its corre- 
sponding anti-Stokes feature (peak la') is observed. At 
around 1 meV a group of three features appears, peak ia 
at 0.8 meV, peak lb at 1 meV, and peak ib at 1.36 meV. 
These transitions were already observed in the previous 
low-energy INS experiment^ and interpreted as follows: 
Peaks Ia and Ib correspond to cold magnetic transitions 
from the 5 = ground state to the first excited 5 = 1 
multiplet, which is zero- field split by magnetic anisotropy 
into its components M — ±1 and M = 0, and peaks ia 
and ib were identified as hot magnetic transitions from 
this first excited 5 = 1 multiplet to the next-higher lying 
5 = 2 multiplet (see inset to Fig. [3]). In addition to these 
transitions, a further peak II at 3.0 meV is observed. 

Figure @|a) presents the INS spectra recorded with an 
incoming wavelength A = 3.2 A on the spectrometer FO- 
CUS. In the 1.5 K data a prominent feature at 3 meV 
is observed on the neutron-energy loss side, which ob- 
viously corresponds to peak II found before in the IN5 
data (Fig. Also, a shoulder near the elastic line at 
ca. 1 meV is observed, which obviously corresponds to 
peak Ib (peak ib is not detected here because of its weak 
intensity and the lower experimental resolution in the 
A = 3.2 A experiment). Two further features are ob- 
served, a very weak feature at ca. 2 meV and a weak 
broad feature at ca. 4.5 meV. The spectrum recorded 
at 75 K shows, apart from a strongly increased back- 
ground due to the excited lattice, no clear features. It 
is therefore reasonable to assume that at this tempera- 
ture predominantly the lattice excitations are observed, 
and the magnetic scattering intensity is distributed over 
all energies. The lattice contribution (on the neutron- 
energy loss side) at low temperatures may thus be es- 
timated by scaling the 75 K data with the Bose fac- 
tor [1 — exp(— E/kBT)]~ 1 , which determines the tem- 
perature dependence of phononic scattering. The esti- 
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FIG. 4: (Color online) (a) Experimental INS spectra recorded 
on FOCUS with incoming neutron wavelength A = 3.2 A at 
1.5 (blue open circles) and 75 K (red open squares). The 75 K 
data were Bose-scaled to yield the estimated lattice scattering 
at 1.5 K (gray solid squares) and were subtracted from the 
1.5 K data yielding the Bose-corrected data (black solid cir- 
cles), (b) S(Q,u) plot of the (not Bose-corrected) 1.5 K data 
shown in panel (a). Intensity is color-coded from blue (low in- 
tensity) to red (high intensity). The labels indicate observed 
transitions, and the asterisks spurion features as discussed 
in the text. Positive energy transfer corresponds to neutron 
energy loss. 



mated lattice contribution is then subtracted from the 
low-tempcrature INS data, a procedure we call Bose cor- 
rection. This approach was used with considerable suc- 
cess in the past i 18 i 34 ' 35 In the Bose-corrected 1.5 K data, 
shown also in Fig.@|a), peak II remains strong, providing 
a strong hint that it is of magnetic origin. 

Also the dependence of the INS intensity on momen- 
tum transfer Q could be studied, which often allows for 
an unambiguous conclusion as regards the origin of INS 
features. The S(Q, io) plot of the 1.5 K data is shown in 
Fig.EJh). Peaks Ib and II are clearly identified, which is 
impressive considering that a non-deuterated molecular 
powder sample was measured at high energies. Magnetic 
and phononic excitations may be clearly differentiated by 
their Q dependence, since for the former the intensity is 
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either strongest at low Q values or typically is maximal 
at around 1.2 A" 1 ^ while for the latter an increase of 
the intensity with Q 2 is expected.— Clearly, the intensity 
of both peaks lb and II is negligible at low Q and passes 
through a maximum between 1.1 and 1.5 A -1 . At higher 
Q values the observed intensity is almost constant. Such 
a Q dependence is characteristic for the magnetic exci- 
tations in AFM molecular wheels Hence, on 
the basis of the temperature dependence, the Bose cor- 
rection, and the Q dependence peaks lb and II are clearly 
magnetic. In contrast, the features at 2 meV and 4.5 meV 
do exhibit their strongest intensity at large Q values, and 
are hence safely assigned to spurious and/or lattice con- 
tributions, which for the 2 meV feature is also confirmed 
by its absence in the A = 4.2 A data, Fig. [3] 

The INS spectra recorded on FOCUS with incident 
wavelength A = 2.26 A at 1.5 and 75 K are displayed in 
Fig. [5fa). In the 1.5 K spectrum peak II at ca. 3 meV 
is again observed. Furthermore, prominent features are 
also observed at around 8.5 meV (peak III) and 12 meV 
(peak IV), and at ca. 6 meV a weak feature is found. 
Peak III is significantly broader than the experimental 
resolution, and appears to consist of two features. The 
75 K spectrum exhibits an increased intensity reflecting 
the strong phononic scattering intensity, and no detailed 
features are observed at the lower energies. However, 
two features are present at energies roughly correspond- 
ing to those of peaks III and IV. A Bose correction as 
discussed before yielded the Bose-corrected 1.5 K data 
shown in Fig.[5ja). Both peaks III and IV are pronounced 
in the Bose-corrected data, which strongly hints towards 
a magnetic origin of this scattering intensity. The right 
shoulder at ca. 9 meV on the high-energy side of feature 

III in the (original) 1.5 K data is significantly reduced 
by the Bose-correction, which suggests a lattice origin of 
this scattering intensity. The S(Q,oj) plot of the 1.5 K 
spectrum is given in Fig. Ob). The peaks II, III, and 

IV are clearly present and exhibit strong scattering in- 
tensity at low Q values, with indications of a maximum 
below 1.5 A -1 , which unambiguously demonstrates their 
magnetic origin. At all energies a significant phonon con- 
tribution is observed at the largest Q values. To gain 
further insight into the nature of the weak feature at ca. 
6 meV, the INS data were analyzed via Q slices, as shown 
in Fig.[5jc). The 6 meV feature is the better detected in 
the Q slices the higher the Q values are, and appears to 
be essentially absent in the Q slice with the lowest Q val- 
ues. It is hence assigned to a spurion or lattice feature. In 
contrast, the peaks II, III, and IV are more pronounced 
in the lower Q slices, again demonstrating their magnetic 
origin. 

To summarize the INS findings, besides the four low- 
energy features la, lb, ia, and ib, which were already ob- 
served before in previous INS experiments^ three fur- 
ther cold magnetic transitions II, III, and IV were ob- 
served in the high-energy regime up to 14 meV. An anal- 
ysis of these transitions using Gaussian fits with sloped 
background yielded the energy positions as peak II: 



(a) 



2000 



1600 



75 K 
1.5 K 

lattice @ 1.5 K 
1 .5 K corrected! 





5 5 

Energy (meV) 



(b) 



> 



5 




0.7 1.0 



4.5 4.9 



(c) 



3* 
c 

CD 



10 

9 
8 
7 
6 
5 
4 
3 
2 
1 




3.7 A < Q < 4.7 A" 
2.7 A 1 < Q < 3.7 A" 
1.7 A" 1 <Q<2.7 A 
0.7 A 1 < Q < 1 .7 A 




4 6 8 10 

Energy (meV) 



12 



FIG. 5: (Color online) (a) Experimental INS spectra recorded 
on FOCUS with incoming neutron wavelength A = 2.26 A at 
1.5 (blue open circles) and 75 K (red open squares). The 75 K 
data were Bose-scaled to yield the estimated lattice scattering 
at 1.5 K (gray solid squares) and were subtracted from the 
1.5 K data yielding the Bose-corrected data (black solid cir- 
cles), (b) S(Q,u) plot of the (not Bose-corrected) 1.5 K data 
shown in panel (a). Intensity is color-coded from blue (low in- 
tensity) to red (high intensity), (c) Q-sliced 1.5 K data, with 
the Q slices as indicated. The labels indicate observed tran- 
sitions and the asterisk a spurion feature as discussed in the 
text. Positive energy transfer corresponds to neutron energy 
loss. 
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3.0(1) meV, peak III: 8.5(2) meV, peak IV: 12.0(2) meV. 



III. ANALYSIS 

A. Spin model for Feis 

The generic Hamiltonian for the description of the 
magnetic properties of a molecular wheel of N = 18 
Fe m spins is a Heisenberg Hamiltonian with nearest- 
neighbor interactions plus a term reflecting the single-ion 
anisotropy of the Fe m ions* 1 ^ In principle, this would lead 
to at least 36 unknown parameters. However, Feis pos- 
sesses a crystallographic C& symmetry axis perpendicu- 
lar to the wheel plane, which leads to a repeating unit of 
three Fe m ions, as indicated in Fig. [IJb). Therefore, the 
Heisenberg part of the Hamiltonian can be formulated 
with at most three different exchange constants: 

L=6 

Hh = ^ Jl&l,l ■ S(,2 + ^Si,2 • Si ; 3 + J^S^ • Sj+i,! , (1) 
1=1 

where L = N/3 is the number of unit cells and I enumer- 
ates the repeating units and has to be understood modulo 
L [for the enumeration of the individual spin sites see also 
Fig. HP>)]. 

The magnetic susceptibility, Fig. [2j is compatible with 
a non-magnetic, total spin S — ground state. In addi- 
tion, the high temperature behavior of x requires that the 
sum of all exchange interactions is antiferromagnetic. 42 
Considering the susceptibility function further, one can 
even conclude that all couplings along the ring have to 
be antiferromagnetic because if any of the exchange in- 
teractions were ferromagnetic, either the ground state 
would possess a total spin of S = 15 or the susceptibil- 
ity would rise to much larger values at low temperatures 
than observed. Hence, we can safely restrict our analysis 
to AFM exchange interactions, J%, J 2 , J3 > 0. A closer 
look at the structure of the molecule suggests the sim- 
pler model where J 2 — J3 — J and J\ = J', since the 
exchange bridges connecting the centers Si t i to Si ,2 and 
Si t 2 to S/,3 are chemically identical and structurally very 
similar, while the bridge connecting centers Si^ to 
is chemically very different (see Fig.[l|. In the following 
we will denote the uniform ring with J\ = J 2 = J3 = J 
as a one- J model, the situation with J2 — J3 — J and 
Ji = J' as a two- J model, and the general case of three 
different exchange constants as a three- J model. 

The interaction with an applied magnetic field B is 
described by the Zeeman Hamiltonian 

N=18 
3=1 

where the g factor is close to 2 for Fe m ions, j enumerates 
the individual spin sites, and /is is the Bohr magneton. 



Previous results suggest that the single-ion anisotropy of 
the Fe m ions, which is modeled by the term 

itf=18 

H D = DJ2(S]) 2 , (3) 

3=1 

is relatively small in magnitude yet strongly affects the 
spin dynamics at very low energies J^- However, it is ex- 
pected to have negligible effect on the excitations at 
higher energies^ Therefore this term is not included in 
our analysis, but this assumption is carefully checked and 
confirmed in Sec. IIII CI 

In order to determine the exchange parameters we ap- 
proach the relevant energy spectrum in three steps: In a 
first step we analyze the spin-wave excitations and com- 
pare them with the observed INS excitations in order to 
narrow down the range of possible parameter values and 
to qualitatively understand the character of the excita- 
tions. This approach is motivated by the L&i?-band pic- 
ture, which connects the higher-energy low-temperature 
excitations to discrete spin-waves, and spin-wave theory 
was indeed able to reproduce the excitations in the AFM 
wheel CsFes with semi-quantitative accuracy^ In a sec- 
ond step we perform large-scale DDMRG calculations for 
quite a number of parameter sets that yield model param- 
eters of high accuracy. Finally, these are compared with 
QMC calculations of the magnetic susceptibility. 

B. Spin- wave calculations 

The one-magnon excitations of Hamiltonian Eq. ([1]) 
can be obtained in the framework of standard spin- 
wave theory (SWT) M- Some modification of the the- 
ory becomes inevitable, however, when studying physi- 
cal quantities such as the dynamical correlation function 
in Eq. (|A3|) . which exhibits divergencies caused by the 
Goldstone modes. The gapped structure of these modes 
in finite Heisenberg clusters and AFM wheels in particu- 
lar can be handled within the framework of SWT by in- 
troducing chemical potentials for the spin sites, yielding 
the so-called modified SWTsi^— We will come back to 
this point in Sec. IIV1 Here we calculate the one-magnon 
spectrum of Eq. (fl|) in a linear SWT approximation; for 
details we refer to Appendix [Cl 

The spin-wave excitation energies of the Feis system 
are displayed in Fig. |6](a) as a function of a shift quan- 
tum number k, which arises from Fourier transforming 
the periodic ring of L unit cells (physically it would cor- 
respond to a wavevector only in the infinite chain). For 
comparison, as guides to the eye, the dispersions of the 
spin-wave branches of a large (L — 300) ring are also 
shown. Three branches wi(fc), u)2(k), W3(/c) are observed 
as expected from the three centers in the unit cell. For 
Feis the k values are restricted to the discrete values 
k = 0, ±;j7r, ±§tt, ^ 7r i and a discrete excitation spectrum 
with 17 energy states is obtained [wi(0) is the ground 
state]. The excitations fall into four energy levels E\, 
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FIG. 6: (Color online) (a) Spin-wave excitation spectrum of 
Hamiltonian Eq. (fTJ) with L — 6 (symbols) and L — 300 
(solid curves) as obtained from SWT for the parameter bo = 
0.074. The energy levels Ex to E 4 are defined by Eq. (fCT5]l . 
The dashed-dotted curves display the spectrum for the one- J 
model, (b) Dependence of Eq. (Q for Ji — J2 = J and J3 = 
J' on the scaling variable y. The two solutions j/i = 0.276 
and 2/2 = 0.687 for &o = 0.074 corresponding to J' /J « 3 
and J'/J « 0.3, respectively, are indicated. Both solutions 
produce identical one-magnon spectra. 



E 2 , E3, and E4 (each is four- fold degenerate because of 
the symmetries in the SWT approximation). The level 
L is formed by the Goldstone mode oj\{k — tt/3), which 
has zero energy in the standard SWTs. It acquires a 
gap however when modified SWTs are used, which in 
the quantum spectrum relates to the gap between the 
5 = ground state and first excited S = 1 multiplet 
(singlet-triplet gap). Within the LlkE-band concept the 
energies in the levels E\ to E4 correspond to the E band, 
and the ground state and k = tt/3 mode to the L band. 16 
The dispersions for the uniform ring or one- J model are 
also displayed in Fig.^a) for comparison (J was chosen 
here such that the maximal energy coincides). 

The SWT calculations allow some useful insight. A 
modulation of the exchange constants away from uni- 
form opens a gap between the branches u>i(k) and LJ2(k), 



which in relative terms shifts E\ down and £2 up. As 
shown in Appendix O SWT implies three relations be- 
tween the reduced energies E1/E3, E2/E3 and E4/E3, 
which depend only on the dimensionless parameter 6q, 



3^1 J2 J3 



(J1J2 + J2J3 + hAY 



(4) 



Detailed inspection shows that the difference between the 
energy levels E3 and E4 does of course vary with varying 
exchange constants, but always remains relatively small, 
and it turns out that it is always significantly too small to 
account for the experimentally observed energy difference 
of peaks III and IV. Finally, the relation E\+E 2 = E4 fol- 
lows from the first two equations in Eq. (IC18|) . The con- 
clusion from all this is that even in the most general case 
of three different exchange couplings, the possible varia- 
tion in the relative energies of the levels is restricted. In 
particular, the three- J model does not introduce signifi- 
cantly more flexibility in adjusting the energies as com- 
pared to the two- J model, despite the additional free 
parameter. Intuitively this is reasonable, since the main 
characteristics which govern the dispersions in Fig. Eta) 
is the maximal energy (or E4) and the size of the gap 
between the energies E\ and E2. This strongly suggests 
that the possible small differences between the exchange 
constants J\ and J2 cannot be resolved in the present 
experiments, and that therefore the two-J model is ap- 
propriate for Fei8- 

Comparing the spin-wave spectrum with the experi- 
mental INS peaks, it is obvious that the level E\ cor- 
responds to peak II (peaks la and lb are related to the 
singlet-triplet gap or level L , which in the present SWT 
calculation is obtained at zero energy). However, at 
higher energies SWT predicts three further levels while 
only two INS peaks III and IV are observed. This results 
in several possibilities for the assignment, which could 
all be ruled out because of discrepancies with experiment 
except one, namely the assignment that E2 corresponds 
to peak III and E4 to peak IV. According to SWT an 
excitation corresponding to level E3 between peaks III 
and IV is expected but not observed. This problem will 
be resolved by the DDMRG calculations discussed in the 
next subsection, which additionally provides scattering 
intensities. 

Due to the restrictions for the energies mentioned be- 
fore (cf. Appendix [Cj , a perfect mapping of the spin- 
wave spectrum onto the positions of the above three 
INS peaks is not possible even in the general case of 
three exchange parameters. A good approximation to 
the INS peaks can however be obtained with the two-J 
model, which is also compatible with the chemical struc- 
ture of Feig. A least-squares fit resulted in the opti- 
mal parameter bo — 0.074, which yields the energy levels 
Ei = 3.17 meV, E 2 = 8.67 meV, E 3 = 10.61 meV, and 
E4 = 11.84 meV, which for E\, E2, and E4 can be com- 
pared to the INS energies peak II = 3.0 meV, peak III 
= 8.5 meV, and peak IV = 12.0 meV. The good agree- 
ment suggests that the exchange constants obtained from 
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this analysis should provide excellent starting values for 
a more refined analysis. 

In fact, from Eq. Q, for any fixed value of 60 there 
are two different pairs of exchange constants which pro- 
duce identical one-magnon spectra (in the two- J model). 
The dependence o f bo on the scaling variable y = 
J 2 1 \J J\ + J 2 + J 3 is presented in Fig. E^b), and the 
graphical solution of Eq. (HJ for the optimal parameter 
bo = 0.074 is also shown. The two solutions correspond 
to J' / J 3 and J' /J ss 0.3, which in the following will 
be referred to as the J' > J and J' < J scenario, respec- 
tively. DDMRG considered in the next subsection allows 
us to disentangle these two cases. 

It is finally mentioned that a good approximation to 
the experimental spectrum cannot be achieved by a uni- 
form ring or the one- J model, see also Fig. Ufa). The 
two-J model is thus both the minimal and appropriate 
model for Feig. This is in accordance with similar find- 
ings for another, but structurally similar Feis spin ring. 4 



C. DDMRG analysis of the experimental data 

In the first part of this section, the magnetic suscep- 
tibility and high-energy INS (peaks II-IV) data are an- 
alyzed. The single-ion anisotropy term is neglected and 
only the Heisenberg Hamiltonian Eq. ([1]) is used. In the 
second part, using the low-energy INS data (peaks la 
and lb) the magnitude and influence of the single-ion 
anisotropy term Eq. (|3]) are investigated. 

Using DDMRG^Sr— (see Appendices \M\ and EJ) the 
INS intensities where first calculated for the one- J model. 
The calculated peaks correspond to transitions from the 
S = ground state to excited states with S = 1, and oc- 
cur approximately at 0.27J, 2.0J, 3.6J, 4.8J, and 5.5J. 
Therefore, as already deduced from SWT, it is not pos- 
sible to find a single J for which more than two INS 
peaks can be reproduced. Furthermore, simulations of 
the magnetic susceptibility (see Fig. [2]) are in marked 
discrepancy with experiment. The one- J model is hence 
disregarded for interpreting the magnetism in Feis. This 
model will however be helpful for the discussion of the 
general physics in AFM molecular wheels in Sec. IIV1 

Simulations of the magnetic susceptibility for the two- 
J model (Fig. [5J showed that for both the J' > J and 
J' < J scenarios the experimental data can be excel- 
lently reproduced. The simulations are in fact not very 
sensitive to small variations of the exchange constants, as 
long as their average remains constant^ Using the infor- 
mation from the SWT calculations and magnetic suscep- 
tibility simulations extensive DDMRG calculations for 
several (J, J') parameter sets were performed that aimed 
at fitting the high-energy INS data. The broadening 
7] = 0.5 meV (cf. Appendix IA 2\ was chosen such that in 
the DDMRG calculation the FWHM of the Lorentzian 
peaks approximately corresponds to the experimental 
resolution (1.1 meV for the A = 2.26 A INS measure- 
ments). For small parameter variations (and for not too 
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FIG. 7: (Color online) INS spectra as obtained by DDMRG 
calculations (solid lines) and comparison to the experimental 
Bose-corrected data (solid and open circles) deduced from the 
A = 2.26 A, 1.5 K INS run. (a) Best-fit simulation for the 
J' > J scenario, with exchange parameters J = 1.42 meV, 
J' = 4.57 meV. The arrows indicate features in the simulated 
spectrum discussed in the text, (b) Best-fit simulation for the 
J' < J scenario, with exchange parameters J — 2.88 meV, 
J' = 1.02 meV. The inset shows simulated INS spectra for 
three slightly different exchange constants used to infer the 
errors in the determined exchange constants. 

small a difference between J and J') the following qual- 
itative results were obtained [cf. Fig. EKb)], which for 
brevity will be discussed for the J' < J scenario (for the 
J' > J scenario interchange J and J'): 

(1) A variation of J' while keeping J constant affects 
the complete spectrum, but the peak positions depend 
approximately linearly on J'. 

(2) Variation of J with constant J' mainly shifts the 
peaks above 6 meV without changing the relative posi- 
tions of those peaks, i.e., the "bandwidth" of the high- 
energy part above 6 meV is not significantly changed. 
The positions and heights of the peaks below 6 meV are 
almost unaffected by variation of J. 

These findings are in accordance with the SWT re- 
sults. Using these trends for the J' > J scenario, the 



best agreement with experiment was obtained for J = 
1.42 meV and J' = 4.57 meV [Fig.[7Ja)]. Five prominent 
peaks are observed in the DDMRG spectrum at the ener- 
gies of ca. 0.43, 3.0, 9.4, 10.9, and 12.1 meV (which can 
in fact contain some nearly degenerate transitions, which 
are not resolved due to the finite broadening r}). The cal- 
culated spectrum is compared with the Bose-corrected 
INS spectrum recorded at a wavelength of 2.26 A in 
Fig. [TJa). The agreement with the experimental peak 
III could probably be further optimized by fine-tuning 
of J and J'. However, in the DDMRG spectra an addi- 
tional peak at about 11 meV is observed, marked by an 
arrow in Fig. [Zta), which obviously corresponds to the 
spin- wave level E 3 discussed in the preceding subsection. 
It is, however, not seen in the experiment, which disfa- 
vors the J' > J scenario as a model for the magnetism 
in Fei8- 

Further efforts therefore concentrated on the J' < 
J scenario. Here also five peaks are observed in 
the DDMRG spectra at ca. 0.43, 3.0, 8.9, 10.1, 
and 11.6 meV. However, the intensity of the peak at 
10.1 meV, which relates to the spin- wave level E$, is rela- 
tively weak. The best agreement of the DDMRG spectra 
with the experimental high-energy data was obtained for 
J = 2.88 meV and J' — 1.02 meV. The agreement is in 
fact very good, see Fig. [7£b), and in particular consid- 
erably better than for the J' > J scenario [Fig. [T^a)]. 
The errors of the determined exchange constants were 
estimated to ~ 0.15 meV for J and ~ 0.05 meV for J'. 
These estimates are based on the variation of the peak 
positions for varied parameters [cf. inset to Fig. Efb)]. 
The smaller error for J' stems from the fact that the ex- 
perimental position of peak II has been measured more 
precisely than the positions of peaks III and IV, and the 
position of the simulated peak II is mainly determined by 
this coupling constant. Within these error bounds, it is 
possible to match the positions of all experimental high- 
energy peaks in the DDMRG simulation. The magnetic 
susceptibility is also reproduced excellently with these 
exchange parameters. 

The consideration of a non-zero but small tempera- 
ture in the simulations would not change the positions 
of the main peaks, but would change their heights a bit. 
Additional hot peaks resulting from transitions from ex- 
cited states would also appear, but because of the singlet- 
triplet gap of about 0.3 meV, which should be compared 
to the temperature of 1.5 K at which the high-energy 
INS data were obtained, the ground state population is 
estimated to be larger than 90 %. Hence, transitions 
from the ground state clearly dominate, and comparing 
the zero-temperature DDMRG and 1.5 K experimental 
spectra is justified. 

Interestingly, the DDMRG spectra for both scenarios 
produce a very weak scattering intensity at ca. 6 meV 
(indicated in Fig. [7] by arrows), where a weak feature is 
indeed observed in the experiment. The analysis of the 
experimental data in Sec. iHl suggested that this intensity 
is of non- magnetic origin. In view of the DDMRG results, 
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FIG. 8: Calculated low-energy spectrum for J = 2.88 meV, 
J' — 1.02 meV as obtained using DMRG with multiple target 
states, (a) Energy spectrum for D = as function of total 
spin S. (b) Energy spectrum for D = 0.03 meV as function 
of the magnitude of the magnetic quantum number M, show- 
ing the ZFS of the S — 1 multiplets. For the lowest S — 1 
multiplet the ZFS is determined by D\. The gray box is cen- 
tered at the center of gravity of the higher-lying multiplets 
and has a height corresponding to the experimental resolu- 
tion (430 fieV) of the data shown in Fig. [t] The ZFS of the 
lowest S — 1 multiplet is 0.7 meV and much larger than the 
splitting of the next two, nearly degenerate multiplets, which 
is 0.3 meV, and smaller than the experimental resolution. The 
observed cold INS transitions are indicated by arrows. 



it could be possible that this feature is, at least in parts, 
i.e., at low momentum transfer Q, due to magnetic scat- 
tering. The present experiments, however, cannot resolve 
this issue. 

In a next step the influence and magnitude of the 
single- ion anisotropy term, Eq. ©, are analyzed. Peaks 
la and lb, observed at energies of 0.3 meV and 1.0 meV 
in the A = 4.2 A INS experiment and in a previous 
workjH were assigned to the transitions from the 5 = 
ground state to the first excited 5 = 1 multiplet, which 
exhibits a Zero-Field Splitting (ZFS) due to the single- 
ion anisotropy of the Fe m ions which can be character- 
ized by an anisotropy constant D\ [see Fig. [5Jb)]. The 
relation between D and D\ was estimated in Ref. by 
extrapolating results on smaller wheels, but the accuracy 
remained unclear, and therewith the accuracy of the de- 
duced value of D in Feig. Therefore we used the ALPS 
DMRG code2£ to calculate the lowest-lying energies for 
H = H H + H D . For J = 2.88 meV, J 1 = 1.02 meV, 
and D — 0.030 meV the low-lying energies reproduced 
excellently the transition energies observed for peaks la 
and lb. Based on parameter variations in the calculations 
and the widths of the experimentally observed peaks, we 
estimate the error in D to be ~ 5 fieV. 

The comparison of the energy level structure with and 
without Hd in Fig. [5] shows that the anisotropy- induced 
splitting is largest for the first 5 = 1 multiplet but much 
smaller for the next two (nearly degenerate) 5=1 mul- 
tiplets. This is in agreement with findings on e.g. the 
CsFes wheels This is a clear indication that the effect 
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of the single-ion anisotropy on the higher-lying excita- 
tions studied in this work can be neglected: The posi- 
tion of the center of gravity is nearly unaffected and the 
splitting is smaller than experimental resolution. It is 
mentioned that the center of gravity is not equal to the 
"center of INS intensity" since the INS intensity is gen- 
erally distributed unequally among the excitations. 

To summarize this section, the experimental mag- 
netic susceptibility and INS data for Feis can very well 
be reproduced with the microscopic parameters J — 
2.88(15) meV, J' = 1.02(5) meV, and D = 0.030(5) meV 
[J = 33(2) K, J' = 11.8(6) K, D = 0.35(6) K]. Alterna- 
tive exchange parameters were extensively searched for 
but found to provide inferior agreement with experiment. 



IV. 



DISCUSSION 



According to the L&_E-band picture, the low- 
temperature excitations in even-membered AFM wheels 
fall into two energy regimes with different character. For 
the wheels Cr g and CsFeg this scenario was confirmed 
both experimentally and theoretically in detail^^ but 
for Feis the picture has been incomplete so far. For Feis, 
the L-band excitations at the lowest energies were ex- 
perimentally investigated in detail in the previous study 
Ref. |l4| by means of low-energy INS (peaks la, lb in 
Fig. [3]) and ultra-low-temperature high-field magnetiza- 
tion and magnetic torque measurements. The L-band or 
the higher-lying elementary excitations accessible at low 
temperatures were carefully studied in the present work. 
However, different models were used in the interpreta- 
tion of the low-energy and high-energy data, and their 
relation is considered now. 

The low-energy experimental dataM could be de- 
scribed extremely well by the effective two-sublattice 
Hamiltonian 



H 



AB 



{S A ) 2 + {S%f 



(5) 



where S,4 and represent the total spins of length 
Sa — Sb = Ns/2 on each of the two sublattices A and 
L, and with appropriate effective exchange constant j 
and anisotropy d. This simple two-spin Hamiltonian is 
suggested by the L&L-band picture, and was established 
as an effective low-energy approximation of the uniform 
ring (one- J) model with an additional anisotropy, or 
H um = Hh + H D with J t = J 2 = J 3 = J. 53 The 
magnetic parameters are related through j — a\J and 
d = b±D, where the coefficients a\ and b\ depend strongly 
on N and s. The underlying assumption is that the 
eigenstates of H uni in the low-energy sectors are well ap- 
proximated by "quasi-classical" spin states of the form 
\S A S B SM) (with S A = S B = Ns/2). Within this 
space H um is equivalent to Hab in first order, yield- 
ing ai = af B with af B = i/N and b l = bf B with 
bf B = (2s — l)/(Ns — 1). However, quantum correc- 
tions modify these parameters significantly, which can be 



accounted for by matching the low-energy states to the 
exact energy spectrum, if the latter is available, yield- 
ing values af 71 and b q ™ M For N = 18, s = 5/2 rel- 
evant for Feis, ffli™ = 0.2721 was extracted from QMC 
calculations, 5 ^ and bf 71 ps 0.07 was determined by extrap- 
olating results of rings with lengths of up to N = 10 i 14 i 53 
Using these results J = 1.64 meV and D = 0.026 meV 
was inferred for Feis from the low-energy data in Ref. 0. 
We recalculated the parameters af 71 and bf n as fol- 



lows: 



qui 



was determined from the singlet-triplet gap 



of the Heisenberg part of H um since it can be obtained 
very accurately by DMRG. For the estimation of 6f m 
the spectra of H un i and Hab were compared for dif- 
ferent values of b\ and D using DMRG and exact di- 
agonalization codes of the ALPS package! 26 ' 55 Our re- 
sults are af™ = 0.2683(1) and bf n = 0.063(2). Hence, 
the refined magnetic parameters J = 1.64(3) meV and 
D = 0.029(1) meV are obtained, with which H un i de- 
scribes very well the experimental low-energy excitations 
in Feis- 

The susceptibility and high-energy INS data, however, 
cannot be reproduced by a uniform ring model as demon- 
strated in Sec. IIII CI The two-J model was required. 
From the susceptibility simulations two equally good pa- 
rameter sets were found (one with J' < J and one with 
J' > J). Based on susceptibility it was not possible 
to prefer one set over the other. Also, spin-wave the- 
ory predicted no differences between these two models 
regarding the excitation energies. The simulation of the 
high-energy INS data using DDMRG revealed differences 
between these two models, and favored the J' < J model 

(Fig. Q. 

Also for the three-J model, or H = Hh + Hd with 
general exchange couplings, the effective two-sublattice 
Hamiltonian Hab is obtained in first-order as the effec- 
tive low-energy approximation. The modulations in the 
exchange (and anisotropy) constants along the wheel are 
effectively averaged out in the L-band states^ suggest- 
ing the average exchange constant J = (J\ + J 2 + Jzj/S- 
However, the parameter a\ is strongly modified. Indeed, 
for the J' < J and J' > J models determined in this work 
the average exchange constants are J = 2.26 meV and 
2.47 meV, while for the uniform ring model J = 1.64 meV 
(all three models yield identical low-energy spectra, e.g, 
energies for peaks la and lb). It is interesting to observe 
that the parameter b\ is comparatively less affected; for 
the models here D = 0.030(5) meV is deduced while the 
uniform ring model yielded D = 0.029(1) meV, which 
agree within the errors. 

The energy spectra of the elementary excitations are 
compared for the J' < J, J' > J, and uniform ring 
models in Fig. [51 where only the Heisenberg parts were 
considered (D = 0). Also, the spin- wave energies as pre- 
dicted by our SWT (Appendix [C| are indicated. Ap- 
parently, the lowest-lying triplet, which belongs to the 
L band, has identical energy in all three models (as ne- 
cessitated by experiment). Furthermore, the next- higher 
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FIG. 9: (Color online) Effect of different coupling constants 
in the Feis wheel. The left panel shows the ground state and 
the excitations in the one-magnon energy sector for the three 
models J = 2.88 meV, J' = 1.02 meV (J'/J « 0.3), J = J' 
= 1.64 meV, and J = 1.42 meV, J' = 4.57 meV (J'/J' « 3). 
The bars represent the energies as obtained from DDMRG 
calculations, and the crosses the results of LSWT. The LSWT 
excitation energies were calculated using Eq. (|C15|) (see Ap- 
pendix[C|. The right panel shows a zoom into the high-energy 
region. 



lying triplet state, corresponding to the energy level E\ 
or peak II, which belongs to the E band, is produced at 
nearly identical energies in all models. The main differ- 
ence is in the further higher-lying states of the E band. 
In the one- J model their "bandwidth" of ca. 2.5 meV 
is obtained approximately correctly, but is predicted at 
significantly too low an energy. Hence, the key exper- 
imental signature for modulated exchange couplings in 
Feis is the large energy gap between peaks II and III. 
Notably, for all three models the simple LSWT can re- 
produce the DDMRG excitation energies in the E band 
semi-quantitatively. This strongly suggests that these 
excitations do indeed correspond to discrete AFM spin- 
wave excitations. The lowest-lying singlet-triplet gap or 
L band is not reproduced by the LSWT (i.e. is calculated 
zero energy) by reasons discussed before in Sec. IIII Bl 

Finally, the excitations of the uniform Heisenberg ring 
shall be considered in detail by comparing the results of 
(D)DMRG calculations to those of various SWTs. The 
discussion parallels that in Ref. 18 for the AFM wheel 
CsFe 8 with N — 8, s = 5/2. AFM systems with disor- 
dered ground states, such as the AFM wheels, represent a 
challenge for any SWT since these start by construction 
from an ordered ground state. As a result, the ener- 
gies of the Goldstone modes are obtained as zero even 
in finite spin systems, and divergencies appear e.g. in 
the magnetization (which actually can be used as an in- 
dication of the absence of order in the considered spin 
model 4 ^). Linear and interacting SWT (LSWT and 
ISWT) are typical representatives ! 44 ' 45 These drawbacks 
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FIG. 10: (a) Dispersion of the one-magnon excitations of 
the uniform Heisenberg ring (one-J model) as obtained by 
DDMRG calculations compared with the dispersion relations 
obtained with LSWT and ISWT (see text for details). The 
LMSWT and IMSWT results are not shown as their disper- 
sions are virtually identical to that of LSWT and ISWT, re- 
spectively, except near k — 0, ±7r/3. The k quantum numbers 
have to be understood relative to the ground state quantum 
number, (b) Oscillator strength calculated with DDMRG 
(stars and numbers) compared with the intensity obtained 
from LSWT. 



can be eliminated by introducing chemical potentials for 
each spin center, which result in finite excitation gaps for 
Goldstone modes (or the singlet-triplet gap) and thereby 
remove the divergencies. This branch of SWTs is de- 
noted as modified SWTs," the simplest of which is linear 
modified SWT (LMSWT) ££M but also a number of in- 
teracting variants exist. Here we use full-diagonalization 
interacting modified SWT (IMSWT) 42 A conceptually 
different approach is to introduce Schwinger bosons and 
treat the resulting Hamiltonian at the mean-field level 
(Schwinger-boson mean field theory) , 50 ' 51 which however 
yields exactly the same excitation spectrum as LMSWT 
and is hence not further considered here. For the N = 8, 
s = 5/2 system it was found that a simple correction of 
the LSWT spectrum, called LSWT+A C , gave the best 
agreement with the exact energies^ In this approach, 
the first-order approximation A c = 4 J/N for the singlet- 
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TABLE I: (D)DMRG and SWT results for the ground-state 
energy E g , singlet-triplet gap A, height of the dispersion 
curve, and mean deviation for a N = 18, s = 5/2 uniform 
AFM Heisenberg ring in units of J. The numerical inaccu- 
racy of the DMRG ground state energy is two orders smaller 
than given. 





E g /J 


A/J 


height/ J 


x 2 


(D)DMRG 


-129.703 


0.2683(1) 


5.50(1) 





LSWT 


-128.835 





4.924 


0.0602 


LMSWT 


-129.082 


0.1228 


4.926 


0.0582 


ISWT 


-129.759 





5.288 


0.0246 


IMSWT 


-129.693 


0.1319 


5.289 


0.0246 


LSWT+A C 




0.2222 


5.146 


0.0316 


ISWT+Ac 




0.2222 


5.511 


0.0066 



triplet gap (or a\ = a^ B ) is added to the excitation en- 
ergies. We will additionally consider ISWT+A C , where 
the estimated gap is added to the spectrum of ISWT. 

For these theories the ground-state and one-magnon 
energies can be calculated analytically for AFM 
wheels4i~— The resulting ground-state energies, singlet- 
triplet gaps, and maximal energies (= heights of the dis- 
persion curves) are listed in Table U which also presents 
the corresponding values for the one-J model Hamilto- 
nian obtained from (D)DMRG, see appendix [Bl Further- 
more the mean deviation of SWT and DDMRG excita- 
tion energies or x 2 = J2ki u ( k ) _ u DDMRG (k)} 2 /(±NJs) 
is given. The excitation energies are also displayed in 
Fig. [Ha). 

The trends for N — 18, s = 5/2 are similar to what 
has been observed for N — 8, s — 5/2: 18 For the ground- 
state energy IMSWT provides the best approximation to 
the DMRG result. The singlet-triplet gap cannot be re- 
produced by LSWT and ISWT, in contrast to LMSWT 
and IMSWT. However, their predictions are about a fac- 
tor two too small, and in particular significantly poorer 
than the simple estimate A c = 4J/N. The calculated 
excitations for LSWT and LMSWT are nearly identi- 
cal as seen by similar mean deviations, as well as that 
for ISWT and IMSWT [LMSWT and IMSWT dispersion 
curves were hence not plotted in Fig. HOT a) for clarity] . 
Thus, as a conclusion, similar to the situation for N = 8, 
s = 5/2 wheels, IMSWT provides the best results also for 
N = 18, s = 5/2 wheels as regards the ground-state and 
_B-band excitation energies, and the singlet-triplet gap is 
underestimated by about a factor of 2. The best descrip- 
tion for the excitations is in fact provided by the simple 
LSWT+A C and ISWT+A C approaches, with a significant 
advantage of ISWT+A C (in contrast to the N = 8 ring, 
where LSWT+A C did better). 

Besides the position of the energy levels also the 
oscillator strengths^ \{S,M\S z (k)\S',M')\ 2 {S z (k) = 
J2j ex P( — ikj)Sj) of the transitions from the ground state 
(S = 0) to the excited S' = 1 states were calculated for 
the uniform N = 18, s = 5/2 wheel using DDMRG. 



The oscillator strength is not identical to the INS in- 
tensity for a specific transition, but both are intimately 
related.— The oscillator strengths are independent of 
specific INS-experiment parameters and are thus better 
suited for general discussions. The results are summa- 
rized in Fig. [TOl b). As expected, the oscillator strengths 
increase with increasing k, becoming maximal at the 
zone boundary, reminiscent to the behavior of Goldstone 
modes at the Bragg points k = n in infinite lattices. 
Interestingly, the oscillator strengths are well described, 
at least on the logarithmic scale used in Fig. HUT b) , by 
the predictions of LSWT^ The LSWT result shows a 
divergency at the Bragg point k = 7r, which reflects 
the underlying assumption of a long-range Neel-ordered 
ground state or zero singlet-triplet gap, as discussed be- 
fore several times. The good agreement of the DDMRG 
and LSWT oscillator strength provides a further strong 
indication that the elementary excitations in N = 18, 
s = 5/2 AFM wheels are indeed discrete spin- wave exci- 
tations, and for the validity of the L&zE-b&nd picture. 



V. CONCLUSIONS 

In summary, we report a combined experimental and 
theoretical study of the magnetism in the antiferromag- 
netic molecular ferric wheel Feis- It is demonstrated that 
nowadays advanced experimental and theoretical tools 
such as inelastic neutron scattering, (dynamical) density 
matrix renormalization group techniques, and quantum 
Monte Carlo together with established methods such as 
spin-wave theory permit a detailed characterization of 
large magnetic molecules with huge Hilbert space dimen- 
sions, such as the molecular ferric wheel Feis. The com- 
bined approach allowed us to determine accurate mag- 
netic parameters as well as to test and rationalize effec- 
tive models such as the L&E'-band picture, which are of 
paramount importance for a deeper understanding of the 
physics in AFM wheels and quantum-spin clusters in the 
mesoscopic regime in general. The concrete result for the 
Feis wheel, whose size is situated at the border between 
microscopic and macroscopic, is that the higher-lying el- 
ementary excitations have the character of discrete anti- 
ferromagnetic spin waves. 
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Appendix A: Methods 
1. Experimental Methods 

The Feis material was synthesized according to the 
procedure described in the literature^ 3 - The magnetic 
susceptibility data was recorded on a microcrystalline 
sample using a SQUID magnetometer (Quantum Design) 
in an applied field of 0.5 T. The date were corrected for 
the contribution of the sample holder. The sample was 
prepared by selecting sufficiently many single crystals 
from the mother liquor, squash them and mount them 
on the sample holder, and to cool the sample down as 
quickly as possible. This procedure typically minimizes 
contamination by magnetic impurities. INS spectra were 
recorded at the time-of-flight disc chopper spectrometer 
IN5 at the Institute Laue-Langevin (Grenoble, France) 
with incoming neutron wavelength A = 4.2 A and at the 
direct time-of-flight spectrometer FOCUS at the Paul 
Scherrer Institut (Villigen, Switzerland) with incoming 
neutron wavelengths A = 3.2 and 2.26 A. For each INS 
experiment a fresh non-deuterated powder sample of Feis 
was synthesized, and filled in a double- walled hollow alu- 
minum can. The weight of the samples were approxi- 
mately 0.5 and 1.8 g for the measurements at IN5 and 
FOCUS, respectively. All data were corrected for detec- 
tor efficiency with a measurement of a vanadium stan- 
dard. At FOCUS also empty can measurements were 
done and used for empty can corrections. Experimental 
resolutions at the elastic line were 165 /LteV in the IN5 
experiment, and 430 and 1100 /j,eV in the FOCUS ex- 
periments. If not stated otherwise, spectra were summed 
over all detector banks. Positive energies correspond to 
energy loss of the neutron. 



2. Numerical Methods 

The INS spectra were calculated as follows. The for- 
mula for the differential INS cross section of powder- 
averaged isotropic system o 36 ' 37 i 59 reads 

cPa k' 



cx — e 



-2W 7^2 



F\Q)S"{Q,u), 



(Al) 



-2W 



where Q = k — k' is the momentum transfer, e 
the Debye- Waller factor, and F(Q) is the magnetic form 
factor of Fe 77/ ions. The scattering function S ZZ (Q, w) is 
defined as 



sin(Qi? 



33' , 



QR 



(A2) 



where Spp (w) is at zero temperature given by 
Sfi,(u) = £<0|£>)<n|3?,|0) 5{hw - E n +E ). (A3) 

n 

For the calculation of Spp (u) the dynamical density- 
matrix renormalization group (DDMRG) technique was 



usedi^Sr— Details of the calculation are given in Ap- 
pendix [B] Strictly speaking the molecule Fe 18 is zero- 
dimensional, but effectively forms a one-dimensional 
chain with periodic boundary conditions and is hence 
suited for DDMRG. However, the calculations are nev- 
ertheless very time-consuming due to the applied peri- 
odic boundary conditions. For the calculations on the 
model with different coupling constants up to 600 den- 
sity matrix eigenstates were kept. For the calculations 
on the uniform Hcisenberg model up to 850 density ma- 
trix eigenstates were kept in order to achieve a smaller 
broadening rj. The truncated weight depends on the en- 
ergies hhj as well as 77, and ranged between 10 -7 and 
10~ 4 , which is a very large, but not unusual value for a 
DDMRG calculation, cf. Ref. EE We checked the results 
for a different number of kept density matrix eigenstates 
and no significant changes of peak positions or heights 
were observed. For the direct calculation of the low-lying 
energy levels the ALPS DMRG code2& was used, and up 
to 3000 density matrix eigenstates were kept. 

The magnetization of Feis wa s evaluated by means 
of Quantum Monte Carlo (QMC) calculations^ 2 ^ 
employing again the ALPS code^ Since the under- 
lying Hilbert space is huge, dimension (2s + 1) N = 
101, 559, 956, 668, 416, we used 10,000,000 sweeps for the 
thermalization and 10,000,000,000 sweeps for the evalu- 
ation of the magnetization for every value of the temper- 
ature. 

All calculations were performed on a BULL supercom- 
puter with 128 cores and 386 GB RAM running ScaleMP 
vsmp. 



Appendix B: Dynamical DMRG 

The basic steps of the DDMRG technique are briefly 
discussed. DDMRG^S r 32 ' 61 is an extension of the stan- 
dard DMRG method^ 7 - -29 ! 62 It is a powerful numerical 
method for the calculation of zero temperature dynam- 
ical correlation functions such as (we set h = 1 in this 
section) 



Ga,b{u) 



-(0|At — 
* E 



1 



-B\0), (Bl) 



■ u> + irj — H 

where |0) denotes the ground state with the energy 
Eq. For the comparison to a spectroscopic experimen- 
tal method such as inelastic neutron scattering, the im- 
portant part of this function is the imaginary part (by 
setting A ~ S z and B — Sp one obtains the function 
ImG A , B (w) = Sff(u), cf. Sec.[A2} 

ImG^H = -(A\— L -\B) (B2) 

7T (E q +oj- H) 2 + rf 

= Y,( A \ n )( n \B)5 v (oj + E Q - E n )(B3) 

a 

where S v (x) is the Lorentzian broadened delta function 
with )im6Jx) = 5(x), and \A) = A\0), \B) = B\0). 
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E n denotes the energy eigenvalue belonging to the eigen- 
state \n). A calculation of the matrix elements {A\n) by 
directly calculating the excited states |n) using standard 
DMRG is possible only for low energies since all energy 
eigenstates up to the desired state have to be included 
as target states in forming the reduced density matrix. 
Many target states, however, decrease the accuracy of 
a DMRG calculation^! The basic idea of the DDMRG 
method is a reformulation of this equation using the so- 
called correction vector, which is defined as^L 



|CM) = 



1 



Eq + oj + irj — H 



\B) 



(B4) 



If one splits the correction vector into \C(u))) = \C v (oj)} + 
i\C l {uj)) with 



|C>)> 



(E + oj - H) 2 + if ' 



\B) 



(B5) 



and 



If possible, we have employed and compared all three 
approaches to obtain and validate the results shown in 
Fig. [TU] However, the numbers shown in Fig. [TU]are not 
numerically exact values. We estimate the relative errors 
of the oscillator strengths to be smaller than 10 % in all 
cases. 



Appendix C: Spin- wave theory for the AFM 
Ji — J2 — Ja Heisenberg Ring 

The one-magnon spectrum of the isotropic Heisenberg 
ring Eq. (|TJ) or three- J model was considered in a first- 
order (next to linear) SWT approximation. Since the 
unit cell of the model contains three half-integer spins 
with s = 5/2, the Lieb-Mattis theorem implies that the 
ground state in the case of AFM bonds is a singlet i 65 i 66 

For the SWT calculations it is convenient to introduce 
spherical coordinates Jr, 9, and <j) through the standard 
relations: 



|C») = g f° V M), (B6) 



a direct calculation of Im Ga,b(u) is possible as 

lmG A , B (oj) = -(A\C t (oj)). (B7) 

|C*(w)) is calculated as the solution of a linear equa- 
tion system within the reduced DMRG basis. The tar- 
get states for the reduced density matrix are |0), \B), 
\C 1 (uj)), and \C r (oj)). For A = B a, different calculation 
of YmGA,B(y>) by reformulating (|B5[) as a minimization 
is also possible^ but we found no significant differences 
between the two approaches. For solving the linear equa- 
tion system we use a simple CG algorithm, 6 ? DDMRG is 
a very time-consuming method because the calculation 
has to be repeated for every oj. However, since calcula- 
tions for different oj are independent, this method is easy 
to parallelize. 

To gain information about the one-magnon dispersion 
relation of the uniform chain we have calculated the dy- 
namical correlation function S z (k,oj), which is defined 
as 



S z (k,oj) 



J2 e lkU ~ 
3,3' 



j) S%{u) 



(B8) 



(k = 2nq/N, q = 0, 1, . . . , 17) and can thus be obtained 
by simply Fourier transforming the S zz , (oj) data^i A dif- 
ferent way to calculate S(k, oj) is to set A = B = S z (k) = 
J2j ex T>(—ikj)Sj for the calculation of ImGA^fw)^ 
However, with this procedure only excitations which con- 
tribute sufficiently to the dynamical correlation function 
can be detected. For the transition from the ground state 
to the lowest three 5 = 1 states it is also possible to di- 
rectly calculate the transition matrix elements and the 
oscillator strengths using standard DMRG (cf. Ref. l64h . 



cos <p sin w = x, 



A 

Jr 
■h 
Jr 

Jr 

Jr = \l J( + Ji + J3, 



— — = cos 9 



(CI) 



where < < § and < 9 < n/2 for AFM exchange 
constants. The radial coordinate Jr appears as an overall 
factor in the Hamiltonian, Eq. (fTJ), and sets the energy 
scale. The spin operators were parameterized as follows: 



Sf a = cos$ 



S?„ = sin$ 



Sf a = y/spi, a , 



1 V2 

2 " 2 [Pl - a 

2 2 [Pl > a 



■ sin$ y/Sqi, a , 
cos$VSqi ta , 



(C2) 



where a — 1,2,3 and / = 1,2, ...,L. The angle $ 
alternatively takes the two values and 7r on the lat- 
tice sites along the ring. The coordinate Q7.a and mo- 
mentum pi.a satisfy the usual commutation relations 

[qi, a ,Pi',p] = ibi,i>5 a ,p- 

In terms of the Fourier transforms 
1 L 

(qk,a,Pk : a) = -7= /Xqi,a,Pl,a) exp (-jkl) (C3) 
V L l = i 

the linear spin-wave theory (LSWT) Hamiltonian reads 



H 



LSWT 



= X! (p- fc ' Mfc ' P fe + ^ k ■ Nfc ' ^ ' 

k 

(C4) 
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where E g = -J R (x + y + z)s(s + l)L, q,t = 
(qk,i,Qk,2,Qk,3) an d Pfc = (Pk,i,Pk,2,Pk,3)- The Hermi- 
tian matrices and Nfc read 



M. 



a/3 



^ x + z x zexp(—3ik) 

x x + y y | , (C5) 

y z exp(3ifc) y y + z 

x + z —x — zexp(— 3ik) 
-x x + y -y I • (CO) 

— zexp(3«fc) — y y + z 



The Hamilton equations of motion for the vectors q\- and 
Pfe take the form 



d ^ k t i\/r - 

—— = J R sM k p k , 

at 

dpk ^ 

—j— = ~JRsN k qk, 
at 



(C7) 
(C8) 



which yield three spin-wave branches with dispersions 
u) a (k) defined by the secular cubic equation 



det[M fc .N fc -^p!l 



0, 



(C9) 



Notice that in this expression the first-order correc- 
tion to the one-magnon energies was introduced through 
Oguchi's renormalization factor For the uniform ring 
or one-J model holds R = 1/2 - 1/tt ps 0.1817. In the 
general case, R is expected to be a smooth regular func- 
tion of the parameters 9 and <p. In principle, one can 
calculate this function in the next-order SWT. However, 
the correction is small (few %)and is hence not employed 
here. For the three- J model of Feis, the spin- wave ex- 
citations are distributed in five energy levels defined as 
follows: 



Lo = ^i(O) = u>i(tt/3), 

Ex = wi (±7r/9) =u>i(±27r/9), 

E 2 = o; 2 (±7r/9) =w 2 (±27r/9), 

E 3 = wa(0) = w 3 (0) = w 2 (tt/3) = W3(7r/3), 

Ei = w 3 (±7r/9) = w 3 (±27r/9). (C15) 

The degeneracies are dictated by the spatial Cq symme- 
try plus the underlying bipartite sublattice structure. 

Since the coefficient ao in Eq. (|C12|) vanishes for k = 
and ±7r/3, the roots at these special points of the Bril- 
louin zone can easily be determined to £i = 0, £2,3 = 1- 
In terms of the energy levels one hence finds Lq — and 



where I is the unit (3 x 3) matrix. In diagonal form the 
spin-wave Hamiltonian Hlswt then becomes 



Hlswt = E g + y^u a (fc)rtfc ;Ct , 



(CIO) 



a— 1 k 



where Hfc iQ is the occupation number operator of the 
(fc, a) state and E g is the ground-state energy: 

E g = - J ' + J l + h Ns (s + l) + \ 

a—1 k 

(Cll) 

It is useful to present the cubic equation for u a (k) in 
the following dimensionless form 



f + a 2 £ 2 + a 1 (, + a o = 0, 



(C12) 



with 



no 



Ax 2 y 2 z 2 



(xy + yz + zx) 3 



sin 2 (3/c), 



a x = 1, 
a 2 = -2. 



(C13) 



The spin- wave energies are related to the three roots £k.. 
of Eq. (IC12|) by the expressions 



Ja(k) = \J J\'h + J2J3 + J3J1 S 1H \fik~a 



R 



(C14) 



E 3 = s (l + VJ1J2 + J2J3 + J3J1, 



(C16) 



and the energy E3 plays the role of an overall prefactor 
in the dispersion of the magnon excitations since 



u> a (k) 



(C17) 



Actually, E3 absorbs the explicit dependence on the pa- 
rameters J%, J2 and J3 apart from the dependence of £k.a 
through ao. 

Useful relations between the energy levels £4, E 2 , E3, 
and E4 follow from Viete's formulas at k = tt/9: £k.i + 

£fc,2 + 6c,3 = —0-2 = 2, £fc,l£fc,2 +C*:,2Cfc,3 +^fe,3^fe,l = &1 = 

1, an d £fc, i6c,2£fc,3 = ^0, where b = -a (7r/9). With 
Eqs. (|C15I) one obtains 



2E 2 3 

El 
b El 



El + El 
E\E\ 

rp2 rp2 rp2 



El, 

Tp2 rp2 
&2 &i " 



Tp2 rp2 
~E i E 1 , 



(C18) 



Equations (IC16[) and (|C18[) are very handy to investigate 
the trends in the one-magnon spectrum with varying ex- 
change constants Ji, J 2 and J3, and led to a number 
of general conclusions which are useful for analyzing the 
INS peaks, see Sec. iHTBl 
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